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Abstract. In this paper we derive the density of the first time T that 
a continuous martingale M with non-random quadratic variation {M). := 
Jq h'^(u)du hits a moving boundary / which is twice continuously differentiable, 
and f'/h G £^[0, oo). 

Thus, this work is an extension to case in which M is in fact a one- 
dimensional standard Brownian motion B, as studied in Hernandez-del- Valle 
(2007). 



1. Introduction 

Patie (2004) [Section 3.6.1], for instance, has pointed out the existing relation- 
ship between first passage time problems and 3-dimensional Bessel bridges. In 
Hernandez-del- Valle (2007), the author uses this observation to derive the density 
of the first time T that a one-dimensional, standard Brownian motion B reaches 
the twice continuously differentiable boundary /. 

The idea is to first transform the first-passage time into a Cauchy problem, via 
Girsanov's and Feynman-Kac's theorems [as has been done before, for instance, by: 
Groeneboom (1987), Martin-L6f (1998), Novikov (1981), and Salminen (1988)]. 
Next, after properly manipulating the corresponding parabolic partial differential 
equation, it is possible to obtain Schrodinger's equation for time-dependent linear 
potential which has been solved analytically by Feng (2001). 

In this paper we derive the density ip of the first time that a continuous martingale 
M with non-random quadratic variation {M). := h?{u)du reaches a moving 
bounday /, i.e. 

T:=inf{t >0|Mt = /(t)}, 
provided that the following two conditions hold: 

/eC2[0,oo) and /'//i e €^[0, oo). 
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The idea, once more, is to first derive and next solve a Cauchy problem, which in 
this case, is of the following form: 

- ^(t, o) = ^f3'{t)av{t, a) + lh^(t)^it, a) 

v{s,a) = 1, where 

(3iu) := ^ .>0. 

To understand better the application of this result, let us start with an example: 
Suppose we want to find the density of first time T that a mean-reverting Ornstein- 
Uhlenbeck process X :— {Xt,t > 0} starting at zero, and with unit diffusion 
coefficient will hit a moving boundary /, i.e. 

(2) T:=M{t>0\Xt^ fit)}. 

It is the solution to this type of problems that we tackle in this paper, since equation 
([2]) is equivalent to: 

^ e^dBu = /(t)e* I 



T := inf O > 



where B is a Brownian motion. 

Furthermore, observe that this type of results may in principle be used for the 
pricing of moving barrier options, for instance, Lo and Hui (2006) study the case 
in which the driving process is a constant elasticity of variance (CEV) process. For 
the case in which the barrier is constant and the process is geometric Brownian 
motion, the reader may consult Karatzas (1997). 

The paper is organized as follows: In Section [2] we introduce the problem and 
state our main results, which we will pursue in the remainder of the paper. Next, 
in Section [3] we introduce some auxiliary results, namely: (1) the density ipa of 
the first time Ta that a continuous martingale M starting at zero reaches level 
a > 0, (2) the transition density of M or its associated Green's function, and (3) 
the transition density of M conditioned to reach level a for the first time at time s 
[see for instance: Imhof (1984) or Revuz and Yor (2005)]. 

We devote Section |4] to verify that the dynamics of the conditioned process 
{M\Ta — s) indeed satisfies the transition density derived in Section[3]by obtaining, 
once again, its Green's function. Finally, after deriving the corresponding Cauchy 
problem ((T|) in Section [5l we solve it in Section [6l 

Remark 1.1. For brevity we have not explicitly stated Feynman-Kac's theorem, 
we thus suggest the reader to consult for instance: Theorem 5. 7. 6, and Definition 
5.7.9 in Karatzas and Shreve (1991). Finally, we would like to mention that the 
tools used in addressing the parabolic p.d.e's studied in this paper, are very much 
in the spirit of Takamizawa and Shoji (2003), and Feng (2001). 

2. The problem and main results 



For expository purposes, we present our main results in this section, and devote 
the remainder of the paper to their proof. 
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Problem 2.1. The main motivation of this work is finding the density tpr of the 
first time T , that a one- dimensional, continuous martingale M with non-random 
quadratic variation 

{M)t := / h^{u)du, t > 0, 



reaches the twice continuously differentiable moving boundary f : 

(3) T ■.= mi{t>Q\Mt^ f{t)}, 
given the following two constraints: 

/eC2[0,oo) and /'//i £ €^[0, oo). 

Namely, the continuous martingale M has the following dynamics: 

dMt = h{t)dBt, Mo = 0, i > 

where B is a one- dimensional standard Brownian motion under some measure P. 

Alternatively, from an application of Girsanov's theorem and the forthcoming 
careful construction (in Section of a Markov process Y which starts at level a 
and reaches for the first time zero at time s {i.e. F is a Bessel bridge type process), 
it follows that Problem 1 2 . II mav be restated as: 

Theorem 2.2. Given that: (i) 

ah'^(t) I I 

(4) ipa{t) = ^= ^^ expJ ^ 

I2t: (^J*h^{u)du 



2QK^{u)du 



is the density of the first time that a continuous martingale M , with quadratic 
variation (M). = h^(u)du, reaches the fixed level a G (0, oo). (ii) The process Y 
has the following dynamics: 

(5) dYt = h{t)dWt + h\t) ( I- - ^] \dt, fo - a e (0, oo), t € [0, s\ 

yYt h\u)du j 

(where W is a Wiener process). Then the distribution of T , defined in equation 
(0^, equals 

t-t 

P(T <t) = 



E expj-^ /3'(u)f„du| 



(6) X exp |— /3(s)a + a j f3'{u)du~- J b{u)(3{u)duj ipa{s)ds 

provided that both: f G C^[0, oo) and f'/h G C^[0, oo) and given that [3(u) :— 
b{u)/h^{u), and b{u) f'{u) for u G K+. 

From ([6]) it becomes clear, that our goal is thus to compute the following expected 
value: 



(7) E 



expj-^ l3'{u)Y^du^ 



We shall prove, in Section [U that in fact equation ([7]) has the following solution: 
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Theorem 2.3. Let: (i) f e €^[0,00) and f/h e C'^[0,(X)), (ii) P'{u) > for all 
u> and, (Hi) the process Y is as in equation 0j in Theorem \2.Si Then 



where: 
Given that 

G{t',y-T',z) 



exp 



(3'{u)Yudu 



G{t, a; s, b)db. 



VMr' - t') 



exp 



exp 



'2{t' - t') 



2(t' - t') 

[for < t < T < s and a,b ^ (0, oo)/, and t' — t :'= h'^{u)du, y :^ a + v{t), 
z := b + v{t) are obtained by solving: 

TTt[t)^P'{t) v'it) = -h^{t)Tr{t) 

n^{T) = -/3'(t) w'(t) = h^{T)n{T). 

Then the Green's function G equals: 

>lbn( 



G{t,a;T,z) ^ ex.p\bTr{T) + aTT{t) + - J h'^{u)Tr^{u) \ G{t' ,y;T' , z). 



Thus: 



Corollary 2.4. The density ipx of the stopping time T defined in equation (0) 
equals: 



ipris) 



id ifa is defined in |7p- 



G{0,a;s,b)db 



X exp < —f3{s)a + a / (3'{u)du 



2 Jo 



b{u)P{u)du > fa{s), s > 0. 



3. Auxiliary results 



This section is devoted to the study and derivation of the transition probability 
of a martingale AI conditioned to reach for the first time some fixed level a > at 
precisely time s > 0. The reader should observe that this analysis will lead to the 
construction of Bessel bridge type process. 

We first begin with a couple of well known ancillary results: 

3.1. Density of the first passage time of AI over a fixed boundary a. 

Observe that if: 

(8) Alt / h{u)dBu, Mo = 0, t > 

Jo 

is a P-martingale and we define the stopping time Ta as 

(9) Ta := inf |s > 0\Ms = a| a > 0. 
Then, the density (pa of Ta under P is: 



(10) 



fa{t) = 



ah^{t) 



: exp ■ 



2^ 



2 y^K^{u)du^ 
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for a > and t>0. 
Alternatively, 

3.2. Distribution of Af . From Feynman-Kac's theorem it follows that the tran- 
sition probability of M, equation is 

(11) G{t,x;s,y) = ^=^^^ exp | - — I 
for t > s. 

3.3. Derivation of the transition probability of the conditioned P martin- 
gale Y. [The techniques employed in this section are very much in the spirit of 
Imhof (1984), but may also be found in Revuz and Yor (2005).] Let M = a — F be 
a time changed P-Brownian motion which starts at 0. 

For x,a E M+ wc introduce the transition densities defined over R+ by 

p{to,a;t,y) := ^ {Yt e dy) / dy , 
= G{to,a;t,y) 

(12) p,{to,a;t,y) := ^iYt e dy,T., > t)/dy, 

where G and are as in equations pT|) . and (O respectively. From the reflection 
principle, it follows that 

r (Ft e dy, > t) := P" (F^ e dy) - ^{Yt G dy, < t) 
= r (F e dy) - P'^(F e d{2x - y)) 
= G{to,a; t, y)dy - G{to, a; t, 2x - y)dy; 

hence 

Px{to,a;t,y) = p{tQ,a;t,y) - p{tQ,a;t,2x - y) 
= G{to,a;t,y) - G{to,a;t,2x - y). 

Furthermore, recall from equation (jlOp the density of the stopping time T^, and let 

(13) f{to,a;t,x) F"" {T, e dt) / dt 

h^{t) { (a^xY 1 

— : exp ■ ^ 



^j27T{Jl h-^{u)duf I 2 /^^ K^{u)du J 
If we define 

^■^'"(•) =P''(-|T^ = s), 

then 

P'^--^^(fi^ Gda;„z = l,...,n) = /(^n, s, x) TTp^^^^^^^ 

f{to,a;s,x) fj^ 

n fjtj, Xj] s, x) . .f W 

i=i J[ti-i,Xi-i;s,x) 
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where = to < < • • • < < s- In our case, Yq — a and Y — 0. Thus 

f{to,x;s,0) f-J^ 

^^ A\ TT /(^i; 2;^; s, 0) , , 

(14) = I I 771 7-^Po(U-i,x,_i;U,x,)dx^. 

T=l J[t^-1,X^-1^,S,0) 

What we will show in the next section is that in fact this transition probability 
corresponds to that of the following diffusion 

(15) dYt := h{t)dWt + h\t) I ^ - j! 1 dt Yo^a, te[0,s]. 

\Yt Jt h?{u)du I 



4. Verification of the dynamics of Y 

We have answered the question: What is the transition probability of M, given 
that it hits the level a for the first time at exactly time s. In this section we will 
verify that the following identiy holds (MtjT = s) = a — Yt under P given that Y 
is defined in (fTS]) . 

To do so, let us first find the transition probability of Y using once more a 
p.d.e approach and then verifying that it equals equation p^ . That is, we must 
solve both the backward and forward Kolmogorov equations and then construct the 
corresponding Green's function, i.e. 



[0,r) X (0,cxj) 



{t,s] X (0, oo) 

Although we will only derive the backward equation in this section (the reader may 
consult the derivation of the forward equation at Appendix [T]) the idea is to reduce 
equations ()16|) and ()17|) to the backward and forward heat equation respectively 
through algebraic transformations. 

Proposition 4.1. [The techniques used in the proof of this Proposition, are very 
much in the spirit of Shoji and Takamizawa (2003).] The solution of the backward 
equation in I116\) is given by: 



, 9 , , / 1 X \ du' , . 

\x h^[u)du 1 ax 

1,2/ , X 

-h^ir) i - - }' ] ^(r, y) 

yy h^u)du J dy 

\y j^h^{u)dul 



u{t,x) = ui{t,x)B{t)exY>{A{t)x'^} 
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where 



A{t) = 



h?{u)du 



3/2 



(18) 



{ui)t{t,x) + ]-h\t){ui),,{t,x) + h''{t)-{ui)^{t,x) = 
2 X 



for X e (0, oo) and t G [0, t) . 

Proof. We will first compute the partial derivatives. If we let the subscripts repre- 
sent derivation with respect to the corresponding variables, then 



u'^{t,x) 



Ut{t,x) 



{ui)^B{t)exp{A{t)x'^} +2xA{t)ui ■ B{t) exp {A{t)x'^} 

{ui)xxB{t) exp {A{t)x'^] + AxA{t){ui).^B{t) exp {A{t)x^} 

+2A{t)ui ■ B{t) exp {A{t)x'^ ] 

+Ax'^A^{t)ui ■ B(i)exp{A(t)a;^} 

{ui)tB{t) eyip {A{t)x'^] + ui ■ Btexp {A{t)x'^} 

+x'^AtUi ■ B{t) exp {A{t)x^ } . 



Substituting in (|16p we have 



.^J^u'-^u'-x'A,u' 
ui B 



2 ui ui 
+h^A{t)u' + 2h^x^A^u' 

1 x 



X h'^{u)dh 



X h'^{u)dh 



{Ul)x I 

u 

Ul 

2xAu', 



cancelling u' 



Ml B 



\h^^-^^ + 2h^xA 



("1)2; 



2 Ml Ml 

+h^A{t) + 2h^x^A^ 



1 



X h'^{u)du 



Ml 
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collecting terms 
(19) 



At + 2h^A^ - 2A 



h'^{u)du 



Ul 



2A 



h'^{u)du 



Ul 2 Ul X Ul 

^+3h'A 

£> 

We first observe from the second line in equation that: 

1 



(20) 



2j^h^{u)du 



Next, if we let r := h?[u)du then dr/dt = —h?[t). Hence, it follows from the 
last line in and equation ([201) that 



1 _ ^ 

bH ^ ^ 2j'h^{u)du 
3 1 dr 

LdB = --dT, 

B 2t ' 



I.e.: 



h'^{u)du 



3/2 



where c is some arbitrary constant. The reader may verify that indeed the first 
term in p9|l vanishes. Alternatively, this will lead to: 



{ui)t{t,x) + ]-h^{t)[ui)xx{t,x) + h^{t)-{ui)x{t,x) = 
2 X 



□ 



as claimed. 

Our next goal is to find a general solution to equation (fTS]) in Proposition 14.1 
and hence obtain: 



Proposition 4.2. The transition probability of the process Y , defined in U5\) . 
equals 



(21) 



Git, x; T, y) = ^^w^^——^H(t, x] r, y) 



f{t,a;s,0) 

where f is defined in equation in Section [2| and 
H{t,x;T, y) 



cxp \ 

^2^[^; h^{u)du] \ I '^[!th^{u)du 



[y - xf 



exp ■ 



{y + xf 



2[^lh^{u)du] 
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[for x,y > and < t < t < s]. In particular H{t, x; t, y) — po{t, x; r, y) where po 
is defined in equation in Section\^ 

Proof. To do so we will introduce onee more a time change r' :— h'^[u)du which 
implies that dr' /dt = h^{t). Next, substituting in equation (fT8|) it follows that 

{ui)r'{T',x) + 1-{ui)xx{t' ,x) + -{ui).j,{t',x) = 0. 

2 X 
Finally, set ui = 1/x ■ U2 to obtain the backward heat equation 

(u2)t'{t',x) + ^{u2)xx{t',x) = 0. 

It is clear, see Appendix [7l at this point that the reduced forward Kolmogorov 
equation U3 will satisfy alternatively: 

1 

2^ 

with the constraint that both x and y are positive. Hence 



-iv3)a'icr',y) + 7;Myyicr',y) = 



H{t ,x;a ,y)= cxp <^ } ~ exp > 



v/2^(a' - r') V I 2(a' - r') j ' \ 2{a' - r') 
or with respect to the backward (i, x) and the forward variables (T, y) as: 
H{t,x;T,y) 



^ 2tt{ h^u)du) \ I 2(/, h^iu)du)j [ 2(/^ h^iu)du)jj 

This in turn leads to transition density of the process Y 

~, , Wy \ h'^(u)du\ 

^' ^' = ^' 

(/7a; |_Jt h^[U)du\ 

as claimed. □ 

It follows that equations and (pi]) are equivalent. The reader may find 
equivalent arguments in for instance: Imhof (1984) or Revuz and Yor (2005). 



5. Proof of Theorem | 

5.1. Change of measure. Recall that if M is the solution of the following sto- 
chastic differential equation 

dMs = h{s)dBs + b{s)ds, Mq = x, 

and 

2 

du < +00 t G [0, +CX3). 

Then under the measure P induced by 



f 




Jo 


h{u)_ 
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the process M is a P-martingale. 

Proof of Theorem \2.Si Given the foUowing continuous P-martingale M 



and letting 



Mt / h{u)dB,, Mo = 
Jo 



inf{t > 0|Mt = /(i)} 

inf{t > 0\Mt - [ f'{u)du = /(O)} 
Jo 



= M{t >0\Mt^ a} under P, 
where the Radon-Nikodym derivative equals dF/dP ~ Z and 



exp ■ 



exp • 



exp 



t rl 



dB,, 



1 



Hu) " 2 



/'(«) 



t PI 



fin) 



h'^iu 
P{u)dMu 



dMu - I 



h{u) _ 

fin) 
h{u) 



du 



du 



(3{u)f'{u)du 



where /3{u) = f'{u)/h^{u) for w > 0. Then 



"{T <t) = 



exp <^ - / (3(u)dMu 



P{u)f'{u)du\l^ 



(T<t) 



E 



exp \ -l3{t)Mt + / (3'{u)M^du --J P{u)f'{u)du \ \T<t) 



expl-P{T)MT + l3'{u)Mudu--J^ (3{u)f'{u)du\l(T<t) 



,-/3(s)a-i/o=/3(«)/'(u)du 



E 



Finally, from equation psp in Section [3] we have: 



exp 



ipa{s)ds 



E 



exp<^ / j3'{u)Mudu 



exp 



T ^ s 



{u)du \ E 



exp 



l3'{u)Yudu 



□ 



6. First passage time 
This section is devoted to the computation of the following expectation 

E expj-^ P'{u)Yudu) 
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which is equivalent to solving the following Cauchy problem: 

du 1 du 

- — (t,a) = -p'(t)au{t,a) + -h'{t)-^(t,a) 

.)x (0,00) 

v{s,a) = 1. 

The idea is to solve simultaneously both the backward and forward Kolmogorov 
equations: 

-^(t,a) = -p'{t)av{t^a) + \K\t)^{t,a) 
(22) +h\t)[-- " )§^(t,a), [O,r)x(0,oo) 

in order to obtain its corresponding Green's function. 

Proposition 6.1. The backward and forward equations in \2'/3j) and i23\) satisfy 
the following relationships respectively for < t < t < s and a,b € (0, oo): 



(p{t,a) = (p^{t,a)—^exp{A(t)a'^] 
a 

^(r, b) = i;^ (r, b)bB{T) exp {i(r)fe2 1 



where 



1 



Bit)^c(^J\' 



3/2 



and 



-^(i,a)+/9'(i)V(<.o) = ik"(i)|^(i,o) 
?f-{r,b) + f/lr)bv'lr.b) = ik2(r)^(r. t). 

Proof. The reader may derive it from the techniques used in Scction|4]and Appendix 

m □ 

The backward and forward Schrodinger equations have the following solution: 
Theorem 6.2. For < t < t < s, and a,b ^ (0, oo), and given that 



A/27r(r' - 1') 



jz - yf \ r {z + yf 
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where t' — t' — h^{u)du, y = a + v{t), z = b + v{t) are obtained by solving: 

Mt)=P'{t) v'{t)^-h^it)Tr(t) 

jT^{T) = -P'{t) v'{t) = h^{T)TT{T). 

Then the following Green's function G: 



G{t, a; T, z) = exp |67r(r) + a-K{t) + \ /i^(w)7r^(u)| G{t' , y; t', 



satisfies both the backward and forward Schrddinger's equations for time- dependent 
linear potential: 

-^{t,a)+p'it)a^\t,a) = lh^(t)^{t,a) 
^ir,b)+P'm\r,b) = lh\r)^{r,b). 

Proof Let (f^{t,a) = A(i, a)e''(*)" and tp^iT,b) = A(T,6)e*(^)^ with n{t) and 7r(r) 
being time dependent variables determined later, we have 

iplit,a) = At(t,a)e'^(*)" + a7rt(t)A(t,a)e'^W" 
ifl{t,a) = Aa(f,a)p7(*)» +7r(t)A(i,a)e^(*)« 

^Ut,a) = Xaa{t,a)e^^''>" + 2n{t)Xa{t,a)e^^*^'' + 7T^it)X{t,a)e^^*^^ 



^l{T,h) = A,(r,6)e*(")^ + 67r,(r)A(T,6)e*(^)'' 
i>l{T,h) = A6(r,6)e*(")'' + 7r(r)A(T,6)e*(^)'' 

i>l{T,h) = A66(T,6)e*(^)'' + 2#(T)A6(T,6)e*(^)'' + 7r2(r)A(T,6)e*(^)'' 



which after substitution 

f^ A- nB' (t^\(t. n^ = ^ 



-Xt{t, a) - a-Kt{t)X{t, a) + a(3'{t)X{t, a) = ^h^{t)X^a{t, a) + h^{t)TT{t)Xa{t, a) 



+ ^h''ity{t)Xit,a). 

and 

Xr{T,b)+bTrr{T)X{T,b)+bp'{T)X{T,b) = ^/i^ (r) Afcfc (t, 6) + /i^ (r)^(r) Af, (t, 6) 

+lh\r)n\T)X{r,b). 

If we perform the time and space transformations y = a + v{t) and z = b + v{t) 
where v{t) and v{t) will be determined later, and set X{t,a) = e{t,y) and A(t, 6), 
i.e. 

Xtit,a) = et{t,y) +v'{t)ey{t,y) 
Xa{t,a) = ey{t,y) 
Xaa{t,a) = eyy{t,y), 
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and 



At(t, 6) = er(T, z) +{i'(r)e2(r, z) 
h{r,b) = iz{T,z) 



then 



-et(i, y) + [y- v{t)] [13' {t) - T,t{t)] eit, y) 



+ [v'{t) + h^it)TTit)] ey{t,y) 



and 



1 



er{T,z) + [z-viT)] p'{T)+Trr{T) e{T, z) = ^h'' {t%,{t, z) 



and setting 



we have 



next, let 



thus 



+ [h^{T)niT)-v'{T)]iz{T,z) 

+ lh\T)7r'iT)eir,z) 



7rt{t)=P'{t) v'{t) = -hHt)^{t) 

j,^ (r) = -P' (r) v' (r) = /i^ (r)# (r) 



-et{t,y) = ^h\t)€yyit,y) + ^h\t)7r\t)e{t,y) 

^t{t,z) = ^h'^{T)e,,{T,z) + ^h^{T)TT^{T)e{T,z), 



€{t,y) = <^2(t,y)exp|-i^ /i^(«)7r^(u)du| 



e(r, 2) = xjj^{T, z) exp < - / h? {u)7:^ {u)du 



-'Ptit,y) = -h^{t)<fily{t,y), 



i^l{T,z) = -/i2(r)C(T,^). 



Finally, introduce the following time changes: t' = Jq h?{u)du and r' = h'^{u)du 
to obtain 



rAr',z) = -^i,{T',z). 
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We now start by constructing a function G which satisfies both equations and the 
constraint that a,b G (0, oo), i.e. 



G{t',y;T,z) 



1 



VMr' - t') 



exp 



2(r' - t') 



exp 



'2{t' -V) 



□ 



7. Appendix 

Proposition 7.1. The solution of the forward equation in |J7[ ) is given by: 
v'{t, y) = vi{t, y)B(T) exp {yl(T)t/^} 



where 



A{r) = 



1 



2j'h^{u)du 

-1/2 



h {u)du 



and 
(24) 



(wi)r(r,y) = ^/i^(r)(ui)yy(r,?;) - /i2(T)-(wi)y(r,y) 



1 



1 



h'^{u)d'b 



vi{T,y) 



for y G (0, oo) and r £ (t, s] 



Proof. We will first compute the partial derivatives. If the subscripts represent 
derivation with respect to the corresponding variables, we have 



Vyir, y) 

v'yy{T, y) 



(wi)j,B(T)exp{A(T)j/2} +2yA(r)wi •B(r)exp{A(T)2/2} 
{v^)yyB{T) exp {A^y^} + AyA{T){v^)yB{T) exp {^(t)^^} 
+2A{t)vi ■ B{t) exp {A{T)y^) 
+4y2A2(T)ui • B(t) exp {A{T)y'^] 
v'^{T,y) = (wi)^B(r)cxp{A(r)a;^} +ui •B^exp{A(T)y2} 
.B(r)exp{A(T)?/2} 

Substituting in (fT7|) we have 



-^^w + —V + y ArV 



2 vi " vi 



+h^A{T)v' + 2hYA^v' 



y J'h'^{u)dv 



y h'^{u)dz 



Vl 

2yAv', 
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cancelling v' 



Vi B 



2 vi 



1 



-fi' 



y h'^{u)du 



y J^h'^{u)dh 
1 1 



Mi 

Vl 



2yA, 



collecting terms 



2h^A^ -Ar + 2A-:j 



+h'y 



Vl 



2A + 



h?{u)du 

1 



h?{u)du 



(«l)r , 1.2(^^1)1 



Vl 2 Vl 



1 {Vl)y 



y Vl 



1 1 

y^ h'^{u)du 



Br 

' B 



h'A 



We first observe that: 



Mr) = - 



2j'h?{u)du' 



Furthermore letting: a = h?{u)du, and observing that da/dr = —h?{T) we have 
that 



1 dB 

B~dT 



1 



2j^"h^{u)du 



or 



~ ~2adT 

i-dB = ---da 
B 2a 



which implies that 



B{t) = c 



J h?{u)du 



-1/2 



and 



(^^i)r(T,y) = \h'^{T){vi)yy{T,y)-h?{T)Uvi)y{T,y) 

^ y 



+h\T) 



1 1 

^~ j'h?{u)du 



vi(t, y) 
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as claimed. □ 
Now set vi = y ■ V2, i-e. 

{vi)t = y{v2)T 

{Vl)y = y{V2)y+V2 
{'Vl)yy = y{V2)yy+2{V2)y 

and substituting into 

y{V2)r = h'^ir)-{V2)yy + h^{T){V2)y 

-h\T){v2)y-h'{T)v2 

+ -h {t)V2 - y , V2 

V Jt h''{u)du 

{V2)t = T^h {T){V2)yy - .s,^, , , V2 



or 



Next set 



V2{T,y) = , V3{T,y) 
n^[u)du 



Thus 



reducing the equation to 

(w3)r(T,y) = ]^h\T){v3)yy{T,y) 

Finally, given the following time change a' := h^{u)du and observing that 
da' /dr = /i^(t), substitute in the previous equation to obtain 

{V3)a'{a',y) = ^{v3)yy{(T',y). 
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